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ABSTRACT 

We consider the dynamics of a mechanical failure induced by a shear stress in a strongly magnetized neutron- 
star crust. We show that even if the elastic properties of the crust allow the creation of a shear crack, the 
strongly sheared magnetic field around the crack leads to a back-reaction from the Lorcntz force which 
does not allow large relative displacement of the crack surfaces. Instead, the global evolution of the crack 
proceeds on a slow resistive time scale, and is unable to release any substantial mechanical energy. Our 
calculations demostratc that for som,e magnetic-field configurations, the magnetic forces cause, effectively, a 
plastic deformation of the crust when the resulting elastic shear stress exceeds the critical value for mechanical 
failure. 
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In two common astrophysical circumstances, a crust of 
a neutron star may be placed under considerable stress. 
Firstly, fast-spinning young pulsars, which arc reducing 
their spin frequencies on a short (~ 1000 years) timescale, 
change the shape of their rotational bulges which deforms 
their crusts. Secondly, in magnctars a ^ 10"'^^ G magnetic 
field exerts a strong twisting force on their crusts. In both 
cases the crust yields, and it is reasonable to assume that 
some of the yield proceeds via explosive ruptures which re- 
lease the crustal strain and produce star quakes, in analogy 
with a similar process in the Earth's crust. 

Neutron-star quakes have attracted significant atten- 
tion in the neutron-star astrophysics folklore. Firstly, it 
has been proposed that the starquakes in the spinning- 
down pulsar may help trigger the sudden spin-frequency in- 
creases known as glitches (see, e.g.. Link & Epstein 1996). 
Secondly, it has been argued (Thompson & Duncan 1995, 
hereafter TD95, see also Blaes et al. 1989) that starquakes 
may be responsible for powering flares in Soft-Gamma Re- 
peaters, a sub-class of magnetars. It is mainly this latter 
suggestion that motivates the present study. Briefly, the 
physical picture is as follows: 

Magnetar emission (see Wood & Thompson 2004 for 



a review) is powered by dissipation of a non-potential 
(current-carrying) magnetic field (TD95, Thompson, Lyu- 
tikov, & Kidkarni 2002). The field exerts Lorcntz force on 
the crust, which is balanced by induced elastic stress. For 
strong enough magnetic fields, Lorentz force may induce a 
stress that exceeds the critical stress of the lattice. This 
leads to breaking of the crust and a release of the seismic 
waves. The seismic waves are coupled to the magneto- 
sphere via strong magnetic field and, after converting into 
the Alfven waves, deposit their energy in a fireball above 
the neutron-star surface. 

TD95 did not specify how exactly the cracking would 
proceed; instead, they assumed that the mechanical fail- 
ure would occur within some finite volume and that part 
of the magnetic/elastic energy contained within this vol- 
ume would be converted into seismic waves (see discus- 
sion in their section 2.2). Molecular-dynamical simulations 
(Horowitz & Kadau 2009, hereafter HK09, Chugunov & 
Horowitz 2010, hereafter CHIO) show that indeed the crust 
breaks down suddenly when it is shear-stressed above some 
critical level. 

A viable alternative theory for the origin of the SGR 
flares has also been proposed (Lyutikov 2003, 2006), This 

theory does not invoke the starquakes as the flares' ori- 
gin. Instead, it proposes that slow untwisting of the in- 
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ternal magnetic field leads to gradual twisting of magne- 
tospheric field lines, on time scales much longer than the 
flare duration. In this picture the magnetosphcric field 
evolves slowly, perhaps due to the Hall drift (Goldreich & 
Reisenegger 1992, Pons & Geppert 2007) or plastic flow 
(Jones 2003). Eventually, the magnetosphere reaches a 
dynamical stability threshold due to increasing energy as- 
sociated with current-carrying magnetic field. Then sud- 
den relaxation of the twist outside the star and associated 
dissipation and magnetic topology change lead to flares, 
in analogy with Solar flares and Coronal Mass Ejections 
(CMEs). In this scenario, the flare itself is produced by 
a rapid reconfiguration of the magnetosphcric field, which 
is assisted by a fast reconnection process, much like what 
happens in the better-studied case of a solar flare. Levin & 
van Hoven (2011) have shown that the external magnetic- 
field reconfiguration during a giant flare would excite the 
large-amplitude torsional oscillations of the NS crusts; thus 
in this picture a flare would preceed the starquake and not 
the other way around. 

It is therefore imperative to understand in detail the 
mechanics of a mechanical failure of the neutron-star crust 
under the strong shear stress. In terrestrial experience, the 
sudden mechanical failure in response to the shear stress is 
common and occurs via a macroscopic shear crack, i.e. the 
breakdown of mechanical elasticity within a thin layer of 
material. The crack's rapid progress is easy to understand 
(e.g., Kostrov 1964): once it forms, the stress at the cracks 
edge is strongly enhanced which facilitates the crack's fur- 
ther propagation along its plane. Often, the efficient crack 
formation requires the solid's ability to form a small void 
at the crack's location, in order for the slippage to be un- 
hindered. Jones (2003) has previously argued that since 
in the neutron-star crust the pressure is greater than the 
shear modulus by two orders of magnitude, the conven- 
tional crack that relies on a formation of the void cannot 
occur. However, molecular-dynamical simulations of HK09 
and CHIO demonstrate a sudden drop in the shear stress 
once the shear of the crust coulomb lattice exceeds a crit- 
ical value of about 0.1. Therefore formation of the void is 
not necessary for the efficient slippage. The situation here 
is somewhat analogous to that of the deep-focus earth- 
quakes. There it is known that the void-aided shear cracks 
do not form; instead it is thought that a localized shear- 
stress-induced phase transition reduces the stress and al- 
lows for the rapid slippage of tectonic plates (Bridgeman 
1945). 

Two remarks are in order. Firstly, molecular-dynamic 
simulations of HK09 and CHIO have not shown shear 
cracks. On the contrary, the mechanical failure propagated 



rapidly through the volume of the simulation domain. 
However, these simulations may not be capable of cap- 
turing macroscopic elasto-dynamical effects: the domain- 
size in these simulations is truly microscopic, of order of 
100 inter-atomic spaces and is orders of magnitude smaller 
than e.g. a mean free path of a phonon. By contrast, the 
localization of the shear failure into a plane does not de- 
pend on the microscopic nature of the mechanical failure, 
but is instead a result of the redistribution of the elastic 
shear stresses once a localized failure is initiated. 

Secondly, in the theory of mechanical failure it is con- 
ventional to distinguish between the brittle and ductile 
cracks (see, e.g., Ashby & Sammis 1990). The former 
is thought to be initiated at a single location and pro- 
ceeds rapidly from the start. The latter begins with a se- 
ries of unrelated micro-cracks and initially proceeds slowly. 
The runaway process is initiated when sufficient number 
of micro-cracks merge to form a propagating macro-crack. 
Some materials under shear stress display plastic fiows 
which may become unstable and concentrate into so-called 
shear-bends (see, e.g., Molinari 1997). All of these pro- 
cesses have an important feature in common: the mechan- 
ical failure and the rapid shearing motion that follows is 
concentrated into a narrow layer. This feature plays the 
central role in our paper. 

In this paper we show that the dynamics of a thin shear 
crack would be strongly affected by the magnetic tension 
inside the crust of a strongly magnetized neutron star. In 
particular we demonstrate that the energy release from 
a thin crack would be strongly suppressed, since the mag- 
netic field provides mechanical connection between the two 
slipping sides of the crack and rapidly suppresses the slip- 
pag 60 The energy released from a thin crack falls short 
by several orders of magnitude to power SGR fiares. The 
strong influence of the magnetic field on the rupture dy- 
namics has no counterpart in the Earth crust, where mag- 
netic field is relatively small and its coupling to the me- 
chanical motion is weak due to the crust's small conduc- 
tivity. We show that for some magnetic-field configuration 
the crust response to the magnetically-induced shear stress 
results in what is effectively a plastic deformation rather 
than sudden rupture, and we construct an explicit model 
of such deformation. 

The plan for this paper is as follows. In sections 2 we 
describe our dynamical solutions for the magnetized thin 
crack, and show how magnetic field suppresses seismic en- 



useful terrestrial analogy is the brittle concrete reinforced with 
metal rods against fracture. Magnetic fields play a role of the metal 
rods. 



2 



ergy release. In section 3 we discuss an explicit exam- 
ple with effectively plastic deformation under the action 
of the magnetar-strength magnetic field. In section 4 we 
conclude. 

2. Dynamics of a thin crack 

In order to estimate the maximal amount of energy re- 
leased by a fracture, we assume that it appears instantly 
as a planar slab of infinitesimally small thickness and an 
infinite lateral extent, and that inside this slab the shear 
modulus and viscosity suddenly become zero. We build 
the simplest-possible model which faithfully represents the 
essential physical aspects of the system. The space is as- 
sumed to be filled with the homogeneous material of shear 
modulus /i, and the material is sheared in the x direc- 
tion; the displacement depends on z only and coordinate 
y is redundant (i.e., the problem is 2-d; cf. e.g. Kostrov 
1966). The crack appears dX t — Q along the plane z = 0, 
thus splitting elastic medium into lower and upper half- 
spaces. The vertical magnetic field component Bz is taken 
to be homogeneous, the horizontal magnetic field compo- 
nent Bx{z) is x-independent, and By = 0. The system is 
assumed to be initially in equilibrium, with 



(1) 



where and T™''^ are the elastic and magnetic stress ten- 
sors, respectively. The above condition can be re-written 
as 

(2) 



Bx J + i;tfc=0 



where 



"^ik — '^f^^ik, (.ik — 2 i^i-k + ^k,i) 



(3) 



Bis the magnetic field , Jis the current density, fi is the 
shear modulus (assumed spacially constant), is strain 
tensor, and is the displacement vector. In the incom- 
pressible medium S^kk — 0, so that aik,k = /^A^^. The static 
magnetic field can be expressed in terms of flux function 
Fix) 

= Bzd,P{z). (4) 
The mechanical equilibrium in the x direction becomes 

gp" = (5) 

(primes denote differentiation with respect to z). Thus, 
the initial displacement is given by 



B 

io = —^P + ki + k2'^z, 
47r/x 



where ki and ^2 are constants. 

The quantities T^l and T™'^^ represent the flow of the 
X component of momentum density in the z direction [see, 
e.g.. Landau & Lifshitz 1956], and thus determine the dy- 
namics of the system. 

When the crack appears at t — along the plane z = 0, 
thus splitting elastic medium into lower and upper half- 
spaces, it hinders the momentum flow since = in 
the craclj^ The momenta of opposite sign accumulate at 
the upper and lower sides of the crack, and they move 
in the opposite direction. In the unmagnetized case this 
motion is indefinite, and all of the elastic energy of the 
medium can be released. However, the magnetic tension 
in the z direction results in a rapid change of the magnetic 
xz-stress within the crack and, as we show below, a rapid 
re-establishment of the initial momentum flow through the 
crack. This results in the slipping motion being stopped 
and the associated energy release being suppressed. 

2.1. Motion generated by the crack 

Let ^(z, t) be the a;-displacement from the initial equilib- 
rium, and let bx{z,t) — SB^/B^^ where SB^ is the change 
in Bx compared to the initial value. The variables are cho- 
sen so that C(z,0) = C(z,0) = K{z,Q) = 0. The force 
balance in the x direction and the induction equations be- 
come 



(7) 
(8) 



(6) 



Here rj is the magnetic diffusivity, and ca — B^j \f^^ and 
Cci = \fpTp are the Alfven and the elastic shear velocity, 
respectively, where \x is the shear modulus. Our notation 
is ^ = dildt, C = dildz, etc. 

Eqns ([7][8]) describe the behavior of the elastic-resistive 
medium. Its normal modes are considered in Appendix [A] 
There two types of modes: the slowly-damped mechanical 
shear waves of the elastic medium with magnetic fleld , 
and resistive diffusion-type modes. Next we demonstrate 
that a sudden crack in a magnetized elastic medium excites 
mostly the slow resistive modes, while the amplitude of the 
magnetized elastic modes is very small. 

A convenient way to solve initial-value problem with 
moving boundaries is to use the Laplace transform. We use 
the following notation: the Laplace transform of a function 

^This assumption is most favorable for the energy release. One could 
just as well assume that takes some lower than the initial but 
non-zero value 
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f{t) is given by 



f{p)= / f{t)cxp{-pt). 



(9) 



The behavior of f{t) at late times t S> 1 is determined by 
the behavior of f{p) for small values of p <C 1. 

We now take Laplace transform of the system of homo- 
geneous equations ([t]) and ([s]) and look for solutions of the 
form exp[A(p)z]. We get the following general solution: 

^{z,p) = A{p)exp{-Xiz) + A{p)exp{Xiz) 

+j3{p) exp(-A2z) + B{p) exp(A2z), (10) 



where 



and 



Ai 



ct P 



Cel V V 
^2 = P/Ct 



(11) 

(12) 



The dispersion relation is obtained by using the smallness 
r)p/cly <C 1. This is an excellent approximation (see be- 
low) for the late-time behaviour of the system. The ho- 
mogeneous solution with A2 corresponds to the familiar 
shear wave propagating under the combined action of the 
elastic and magnetic restoring forces, while that with Ai 
corresponds to the "diffusion wave", i.e. the harmonic per- 
turbation damped by magnetic diffusion in the z-direction. 

We now match the general solution above to the appear- 
ance of a crack at z = 0, t = 0. In the upper half-plane 
z > only the outgoing waves exist, so A ~ B = 0. The 
values of the coefficients A{p) and B{p) is obtained from 
two boundary conditions at the crack surface. 

The first boundary condition is straightforward: the 
elastic stress at the boundary is zero, so the shear at the 
boundary has to be zero and therefore 



e'(o+,t) = -eo(^ = o). 



(13) 



Here £,o{z) is displacement of the crust from the position 
of zero elastic stress, before the crack appears. In the case 
that the elastic stress balances out the magnetic one, we 
have 



eo = - 



BzBx 

ATlpL 



Bx 



(14) 



where = \f[ijp is the velocity of a purely elastic shear 
wave. Taking the Laplace transform of Eq. ( 13 1, we obtain 



^"'(0+,p) = -^eo(^ = 0). 



(15) 



From here onwards we shall write the shorthand of for 
C^(z = 0). 

The second boundary condition is obtained by consid- 
ering the singularity at the crack, ^' = 2^(0+)(5(z). Inte- 
grating Eq. (|8| across the boundary, using the fact that 
is finite, and the symmetry h'^(z) = — 6^(— z) we get 



7y6;(0+,t) + e(0+,i) =0 



(16) 



Taking the Laplace transform of this equation, we ob- 
tain 

776L(0+,P)+K(0+,P)=0. (17) 
From Eq. ([7| we get 

6^ = (l/ci)[p2e-4n. (18) 



Substituting Eq. (18 1 into Eq. ( 17 1 we get at the boundary 

z = 0+: 



P 



PV 

'-A 



1 



(19) 



Substituting Eq. (|10| (with A = B ^ 0) into Eq. ([19|) we 
obtain (neglecting rjp/c^^): 



B ^ -fA. 



(20) 



dition in Eq. (13), we get 



Substituting this and Eq. (|10|) into the first boundary con 

1 



where 



Po 



4 4 
4v 



(21) 



(22) 



We consider two useful limiting cases of the above equa- 
tion (the general solution will be obtained later in this sec- 
tion): 

Case (i): For unmagnetized crust, ca = Po = 0, and one 
gets 



B{p) = c,,i'Jp\ 

Mp) = 0. 



(23) 



The inverse Laplace transform of B{p) is given by 

B[t) = etc; (24) 

and the corresponding displacement in the upper half- 
plane z > is given by 



e t 



(25) 
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where 9(i) is the Heaviside function. This solution rep- 
resents the shear wave launched from the suddenly intro- 
duced unmagnetized crack. The velocity jumps discontin- 
uously across the crack, 



e(0+,t)-e(0_,i)-2cei^[„ 



(26) 



see 



Fig.g 







V 










Fig. 1. — Examples of the displacement after the crack 
in a non-magnetized material {Left Panel) and magnetized 
material {Right Panel) for different times < < ^2 < 
^3. Initial profile is — tanhz, both shear and Alfven 
velocities are unity, r] = 0.01. M t = the crack appears 
at the z = plane, launching two shear waves in the non- 
magnetized case and two shear and two resistive waves in 
the magnetized case. In the magnetized case the shear 
waves have very small amplitude, oc and do no carry 
large energy flux. The resistive waves are limited to a 
narrow region Az ^ /. 

If magnetic field is present, however, this fast slippage is 
quickly stopped by the rapidly building magnetic tension, 
and all further slippage is controlled by magnetic-field dif- 
fusion. This leads us to consider another limiting case. 

Case (ii): Of major interest for us is the solution at 
late times, t ^ 1/po- To these, only p <^ pq contribute 
substantially and in this limit we have 



B{p) 



3/2 

; — 2 ^ ' 



and 



A{p) = §B{p) 



Substituting this into Eq. ( 10 ) , we obtain 



^{z,p) = iB{z,p) +£.a{z,p), 



(27) 



(28) 



(29) 



where 



and 



iB{z,p) = ^=^^exp 



U{z,p) = — Co-^exp 

Taking the inverse Laplace transform, we obtain 
2 




(30) 
(31) 



CB{z,t)=C, 



0^2 



c\ct ^/7f 



.jTl-{t-z/ct)Q{t-z/ct), (32) 



dU{z,t) 
dz 



and 

U{z,t)^^'oi 



-Co ■ Erfc 
exp(-[z//]2) 



z/l] ■ Erfc(z//) 



(33) 



(34) 



where 



2cci 
0.02 



rjt 



Ct 



V 



10 ^cm^s 1 



1/2 



t 

lOs 



1/2 



(35) 



(see Fig. [T]). 

The two contributions have distinct physical character. 
The ^B{z,t) represents a shear wave which is launched 
by the slippage at z = and propagates with the speed 
Ct- It is this contribution which is responsible for the re- 
lease of seismic energy into the crust. On the other hand, 
the ^a{z, t) represents the magnetic-diffusion-type process 
which operates locally within distance ~ I given by the 



equation Eq. ( 35 1 from the crack. Within this thin layer 
the elastic stress of the crust is substantially reduced. 
The surface z = 0+ moves according to 



az = o+) = eo^V^ 

Cj^Ct 



Thus, within a very short time. 



t 



2^ c\ 



(36) 



(37) 



the slippage velocity falls below the elastic shear velocity: 
magnetic field effectively stops the relative motion of the 
crack's surfaces. 

For typical neutron-star crust parameters (Chamel & 
Haensel 2008), one evaluates 



Pa 



lO^is-i 



1025s- 



Ct 



lO^cm/s 



(38) 
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where 



/? = <?aIcI 



(39) 



For magnetar-strength fiels, /3 ^ 1 and all times of interest 
satisfy t 3> 1/po- Therefore the approximations for the 
case considered here hold extremely well in a magnetar. 
The energy flux going out in seismic waves is 



dxdydt 



= 3 X 10' 



rfl2 



1014 108 



Ct 

0.1 



t 

V erg 
lO^^i cm^s ' 



(40) 



where all the quantities in the expression above are ex- 
pressed in the cgs units. 

Taking an optimistically large crack area of lO^^cm^, 
we get the typical released energy of ^ lO^^erg for the 
magnetar-strength field with (3^1. This is ^ 12 orders 
of magnitude short compared to typical weak SGR flares 
(Woods & Thompson 2006). Clearly, thin shear cracks 
cannot be responsible for the flares. We finish this section 
with 2 remarks. 



One can evaluate the inverse 



(211, thus finding the motion 



Mathematics remark. 

Laplace transform of Eq. 
due to appearance of the crack without considering limit- 
ing cases. It is easiest to first evaluate "seismic velocity" 
B, whose Laplace transform is pB{p). The latter could be 
calculated by noticing that the inverse Laplace transform 
of 



I{P) = 



1 



(41) 



is given by 



fit) 



, P + l 
e* l-ErffVi 



dk 



(42) 



This can be checked directly from the definition of Laplace 
transform (the integrals over first t and then k are elemen- 
tary). The "seismic velocity" is then give by 



B{t)=cti'j{p^t). 



(43) 



It is straightforward to check that the limits t ^1 /pq (free 
slip) and t » 1 /po (magnetic-diffusion-controlled slip) are 
recoverecQ 



^For X 2> 1, with high degree of precision f{x) = ^/Y/Jjrx). 



Physics remark. We have considered the configuration 
where there is a substantial magnetic field in the direc- 
tion perpendicular to the crack. One could ask whether a 
crack direction could adjust itself in such a way that the 
crack surface would be tangential to the field lines. From 



Eq. ( 40 ) , we see that the Bz component would have to be 
smaller than the average B hy a factor greater than 10^ 
in order for the crack to release sufficient energy for the 
flare (this is because (3 (x Bl). It seems very unlikely to us 
that this type of fine-tuned alignment could be efficient on 
a large scale, unless there was an elasto-dynamical mech- 
anism that would drive the crack into alignment with the 
field. We have so far failed to identify such a mechanism. 

3. Plastic flov^fs? 

We now construct an explicit scenario where changing 
magnetic field configuration induces effectively a gradual 
plastic flow in the constant-density crust. We shall assume, 
as before, B^ — const and consider B^ that is initially 
increasing with time 



B^ 



Ct[cosh(z/L)] 



(44) 



due to some non-MHD process (e.g.. Hall drift) the 
specifics of which are not important for our purpose. Here 
C is some constant and L is the vertical length scale. 
We assume that the xz-component of the magnetic stress 
is balanced precisely by the xz-component of the elastic 
stress, i.e. 



T^mag 
xz 



(^) 



BMz) 
47r 



(45) 



Let T°l-^^ be the critical shear stress at which the lattice 
begins to slip. According to HL09, T°^\t ~ 0.1/x. Once the 
external magnetic field exceeds the critical value at height 



^crit 

47rB, 



(46) 



the crust begins to break and forms a shear crack at that 



location. From Eq. (44 1, this occurs when 



cosh(zcrit/i) 



Ct 



B 



crit ' 



(47) 



The vertical velocity with which the boundary of the crit- 
ical domain moves is given by 



dt ~ 2 \ B"H 



[smh{z/L)]-\ (48) 
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Both the elastic stress and the are at the crack's 
location and are reduced within distance 6z ~ ^/ri5t, where 
5t is the time interval that has passed from the crack's 



formation; see Eq. (35). The time interval 5t between the 



formation of two sequential cracks is given by 
so that 



(49) 



5t ^ rj/vl,,, (50) 
and the separation between the cracks is 

SI ~ T]/?^^!- (51) 

Some numerical estimates are in order. Assuming that the 
magnetar field evolves due to Hall drift on a timescale of 
~ lOOyr, one gets 



Vcrit 



1km 
lOOyr 



3 X 10"^cms~\ 



(52) 



and therefore 



SI 
St 



0.3cm 
lO^s 



day. 



(53) 



It is on the latter timescale that the magnetic field and 
the elastic shear ^' — relax to zero within the stressed 
domain where multiple cracks have appeared. 

4. Discussion 

The neutron star crust cracks differently from that of 
the Earth. In the latter, the quakes are caused by a sud- 
den elastic stress release along a 2-dimensional surface, the 
shear crack, which causes slippage and release of mechan- 
ical energy. The magnitude of the release is much greater 
than the amount of elastic energy stored within the vol- 
ume of the crack. By contrast, in the neutron-star crust 
magnetic tension strongly suppresses the slippage and the 
energy release. Infinitely thin cracks release energy via 
magnetic-field diffusion, too slowly to be able to contribute 
to the energetics of magnetar flares. If SGR flares are as- 
sociated with starquakes which are generated by a crystal 
failure of the stressed crust, these constraints place strict 
requirements on the geometry of the failures. 

We have so far considered simple field geometries, and 
thus have not proved that plastic fiow is the generic macro- 
scopic response of the crust. One can imagine situa- 
tions in which mechanical failure in some region will cause 
magnetic-field reconfiguration that will in turn cause su- 
percritical stress in another part of the crust, thus leading 



to a runaway process. This possibility has to be investi- 
gated further using specific magneto-elastic calculations. 
For now we remark that it is an open question whether 
magnetar flares are associated with the crust failures; as 
was already explained in the introduction, viable alterna- 
tive model is that the flares are caused by reconnection 
events in the magnetosphere, in a manner similar to that 
observed in the solar flares (Lyutikov 2003, 2006). There- 
fore, in our opinion the inferences about flare statistics 
based on phenomenological models for the crust crystal 
failure (e.g., Perna & Pons 2011) have to be made with 
caution. 

We thank Andrei Beloborodov, Andrei Chugunov, 
Charles Horowitz, Rimma Lapovok, Bennett Link, and 
Jay Melosh for stimulating discussions. 
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A. Normal modes in magnetised elastic medium with resistivity 

In this Appendix we discuss the normal modes in the elastic- resistive medium. In Eqns (7][8|, eliminating in favor 
of ^, the equation for the normal modes becomes 

{dt - vd!) {dl - cldl) i = v\dtdl^ (Al) 

This is a linear equation with constant coefiicients and it can generally be solved by separation of variable. Assuming 
C = F{t)G{z), we find 

Fd) f ^F' F"\G" 2 G^-*' 



^^]^-^<i^=^ (A2) 



F ' F J G ' " G 
For harmonic spacial oscillations cx e^'^*+'^^, the dispersion equation is 

cj^ + ik^Tjoj^ - k^cjuj - iclirjk^ = (A3) 
In the limit yy^cs — >■ 0, the normal modes are are 



cjs = -i-^^rj (A4) 

The modes {1,2} are just the resistivity- modified normal modes of the elastic medium with magnetic field . They become 
over-damped (real part smaller than imaginary) for 

k > (A5) 

The mode 3 is the resistive mode. 
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